A transcriptome dataset for gonadectomy-induced changes in rat spinal cord

Circulating sex steroid hormones are critical for neural function and development of neuroplasticity in many regions of the central nervous system. In the spinal cord, our knowledge of steroid hormone influence mostly derives from mechanistic studies of pain processing in dorsal spinal cord circuits; less is known regarding hormonal influence of ventral spinal motor function. Gonadectomy (surgical removal of the testes in males and ovaries in females) rapidly and persistently reduces circulating sex steroids in both females and males, providing a means to interrogate the role of hormones on neural function. Here we provide a next-generation RNA sequencing (RNA-seq) data set to evaluate the impact of gonadectomy on the transcriptome of ventral spinal cord tissue of adult female and male rats.


Background & Summary
Our understanding of how sex steroid hormones influence central nervous system function across the lifespan has rapidly advanced over the past 30 years. The principle steroid hormones estrogen and testosterone are now recognized to exert pleiotropic effects on virtually all areas of the brain and spinal cord 1 . Primarily synthesized in the gonads and adrenal cortex, steroid hormones are released into the blood circulation, readily cross the blood brain barrier, and influence neuronal signaling through nuclear and non-nuclear receptor mechanisms 1-3 . In the spinal cord, steroid hormones directly modulate neurotransmitter release and synaptic neurotransmission 4 , facilitate myelination 5,6 , provide neuroprotection 7-10 , support axonal regeneration 11 following spinal cord injury, and mediate the induction of neuroplasticity 2,12-18 . This knowledge comes primarily from studies of pain processing in the dorsal spinal cord, hippocampal learning and memory, and from rodent models of neuroprotection following spinal cord injury. Less understood is the direct role steroid hormones play in ventral spinal cord function and spinal motor output.
Gonadectomy, the surgical removal of the primary sex organs (ovaries in females; testes in males), induces a rapid and sustained reduction in circulating steroid hormones 2,17,19 . This relatively simple procedure is routinely used in basic science studies to ascertain the role of steroid hormones on function, particularly in the CNS. However, gonadectomy-induced reductions in steroid hormones are likely to alter the fundamental properties of the neural system being studied, complicating comparisons with gonadally-intact groups. This study used next-generation RNA sequencing (RNA-seq) to evaluate the impact of gonadectomy on the transcriptome of ventral spinal cord tissue of adult female and male rats. These data provide valuable, sex-specific gene expression profiles for the ventral cervical spinal cords of female and male rats and suggest unique, sexually dimorphic transcriptomic changes resulting from gonadectomy. geries to reduce levels of circulating sex hormones. This is a common procedure to study the impact of circulating sex hormones on neural function and consists of removing the ovaries bilaterally in female rats (ovariectomy) or removing the bilateral testes (castration) in males. The castrated male rats in this study are referred to as gonadectomized (GDX) and females as ovariectomized (OVX).
At least 2 hours prior to surgery, all rats were provided Buprenorphine SR-LAB (1 mg/kg; s.q. ZooPharm; Fort Collins, CO, USA) for long-lasting pain control. For OVX surgeries, female rats were anesthetized with isoflurane in a closed chamber and then maintained via nose cone (2-3% in 50% O 2 ) for the duration of surgery. Adequacy of anesthesia was confirmed by the lack of response to toe-pinch and eye-blink reflexes. Following removal of fur and skin cleansing with chlorhexidine (Prevantics; PDI Inc. Orangeburg, NY, USA), bilateral incisions were made through the dorsolateral skin and muscle layers to expose the ovarian fat pads. The ovaries were exteriorized and separated from the uterine horns using electrocautery to minimize bleeding. Uterine horns were then replaced into the abdominal cavity. Muscle layers were approximated and sutured with 4-0 absorbable suture (Covidien, Mansfield, MA, USA), and the skin incision was closed with wound clips (9 mm autoclips, MikRon Precision Inc. Gardena, CA, USA). Control rats (n = 6) were ovarian-intact and age-matched to OVX rats.
The GDX procedure in males was performed using similar surgical preparations. Male rats were anesthetized with isoflurane in a closed chamber and then maintained via nose cone (2-3% in 50% O 2 ) for the duration of surgery. Adequacy of anesthesia was confirmed by the lack of response to toe-pinch and eye-blink reflexes. From a supine position, fur was removed over the bilateral testes and skin cleansed with chlorhexidine. Testes were exteriorized via small, bilateral scrotal incisions. The testes were removed using electrocautery to minimize bleeding and vas deferens were replaced back into the scrotal sac. The skin incisions were closed using veterinary surgical glue (Vetbond, 3 M, St. Paul, MN, USA) with special care taken to avoid the penile or anal openings. Control rats (n = 5) had intact testes and were age-matched to GDX rats. Following OVX and GDX surgeries, rats were recovered in temperature-controlled chambers until mobilizing independently, and then were returned to their home cages.
Ventral cervical spinal cord tissue isolation. Rats were deeply anesthetized with isoflurane in a closed container, transcardially perfused with ice-cold saline (pH 7.4), and the spinal cords immediately extracted from C3-C6 (Fig. 1). The dorsal half of the tissue segments were discarded to isolate only the ventral spinal cord. Samples were quickly transferred to ice-cold RNAlater solution (Invitrogen) and refrigerated at 4 °C for 48 hours. Cords were subsequently removed from RNAlater solution, placed into fresh 1.5 mL conical tubes and stored at −80 °C until processed for RNA extraction (Fig. 1). The ventral cervical spinal cord in this region contains upper motor neurons of UE musculature and the phrenic motor pool, a column of motor neurons innervating the diaphragm muscle via the phrenic nerve. Spinal cord tissue from rats receiving castration or ovariectomy was processed 2 weeks post-surgery (Fig. 1), to align with recent studies indicating that respiratory function and the expression of phrenic neuroplasticity are impacted by gonadectomy in rats 2,16,17 .
RNa isolation and sequencing. All samples were sent on dry ice to the University of Minnesota Genomic Center for RNA isolation, library preparation and sequencing. RNA was isolated using the RNeasy mini kit (Qiagen Corp., Santa Clarita, CA) according to manufacturer's instructions. The purity of RNA was checked with Nanodrop (Thermo Fisher Scientific). All RNAs showed a 260:280 ratio of >2.00 (Table 1) RNA isolates were quantified using a fluorimetric RiboGreen assay and RNA integrity was assessed using capillary electrophoresis (Agilent BioAnalyzer 2100) to generate an RNA integrity number (RIN). All samples had RIN values ≥7.8 (avg. RIN 8.2; Table 1) and at least 1 ng total RNA. Sequencing libraries were generated www.nature.com/scientificdata www.nature.com/scientificdata/ using a SMARTer Stranded Total RNA Seq v2 -Pico Mammalian kit (Takara Bio). Briefly, >1 ng of total RNA were fragmented and then reverse transcribed into cDNA using random primers, with a template switching oligo incorporated during cDNA synthesis to allow for full length cDNA synthesis and retain strand specificity. Illumina sequencing adapters and barcodes were then added to the cDNA by PCR, followed by cleavage of ribosomal cDNA. Uncleaved fragments were then enriched by PCR for 12-16 cycles. Final library size distribution www.nature.com/scientificdata www.nature.com/scientificdata/ was again validated using capillary electrophoresis and quantified using fluorimetry (PicoGreen). Indexed libraries were normalized and pooled for sequencing. Libraries were then loaded onto a NextSeq. 550 (single read, Illumina) cartridge, where clustering occurred on-board the instrument. After clustering, sequencing was commenced using Illumina's 2-color SBS chemistry. Following sequencing, Base call (.bcl) files for each cycle of sequencing were generated by Illumina Real Time Analysis software. Primary analysis and de-multiplexing were performed using Illumina bcl2fastq software. The end result of the bcl2fastq workflow was de-multiplexed FASTQ files for subsequent analysis. 151 bp fastq single-end reads (mean read depth ≥40 million per sample) were trimmed using Trimmomatic (version 0.33) enabled with the optional "-q" option; 3 bp sliding-window trimming from 3′ end required minimum Q30. Quality control of raw sequence data for each sample was performed using FastQC (Supplementary Figs. 1-22). Read mapping was performed via Hisat2 (v2.1.0) using the rat genome (Rnor_6.0.98) as reference. Mapping results are provided in Table 2. Average read-pairs per sample was 44.4 M with an average input fragment size of 188 bp ( Table 2). Gene quantification was performed using Subread featureCounts for raw read counts ( Fig. 2A). Sample variance was identified using principle component analysis (Fig. 2B). Clear separation between male and female samples indicate highly sex-specific expression profiles, and samples tended to cluster based on treatment (i.e. gonadally-intact vs. gonads removed; Fig. 2B). The impact of gonadectomy on gene expression is presented in Fig. 3. Differentially expressed genes (DEGs) were identified using the edgeR (negative binomial) feature in CLCGWB (Qiagen, Redwood City, CA) using raw read counts. The frequency distribution of DEGs relative to their fold-change are plotted at three levels of statistical stringency (p < 0.1; p < 0.01; p < 0.0001; Fig. 3). Removal of the gonads appeared to have a more robust effect on the spinal transcriptome of male rats as a total of 1642 genes were significantly up or down regulated at p < 0.1, while 273 genes remained significantly different at p < 0.001 (Fig. 3A). Of these, 191 genes were significantly upregulated and 82 downregulated. Female gene expression in the ventral cervical spinal cord was minimally affected by removal of the gonads. Only 114 genes were significantly changed following OVX at lower statistical stringency (p < 0.1), with 6 DEGs remaining significantly upregulated at p < 0.001 (Fig. 3B).

Data Records
Tables 1, 2 provide specific information related to each sample used in this study. Complete RNA-seq data were deposited in the NCBI's Gene Expression Omnibus (GEO) database (GSE202381) 20 .  Differential Gene Expression following GDX. Removal of the primary sex organs induced transcriptomic changes in the ventral cervical spinal cord. The frequency distribution of DEGs relative to their fold-change is plotted at three levels of statistical stringency (p < 0.1; p < 0.01; p < 0.0001) in male rats following GDX (A) and female rats following OVX (B). Included are the top 10 most differentially expressed genes at p < 0.001 for males, and the 6 genes upregulated in females at p < 0.001.